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Introduccion 



Desde 1935, cuando Einstein, Podolsky y Rosen presentaron un 
trabajo cuestionando la completitud de la mecanica cuantica, se reco- 
noce que el entrelazamiento de sistemas cuanticos es una de las pro- 
piedades mas intrigantes de la mecanica cuantica. En los ultimos 20 
anos se ban desarroUado tecnicas, tanto experimentales como teori- 
cas, que permiten utilizar esta propiedad como recurso para resolver 
problemas computacionales con eficiencias muy superiores a las que 
clasicamente se habian demostrado insuperables. En anos recientes 
se ban logrado tambien avances muy importantes en la comprension 
y clasificacion del entrelazamiento cuantico, aunque subsisten aun 
muchos puntos importantes por resolver e investigar. 

Estudio en este trabajo la capacidad de generar entrelazamien- 
to de una cadena de espines con acoplamiento de Heisenberg XY y 
un campo magnetico uniforme a partir de un estado inicial en el que 
los espines estan completamente alineados. Se encuentra que la ca- 
pacidad de generar estados entrelazados no muestra un comporta- 
miento monotono con el campo presentando, en cambio, "plateaus" y 
resonancias. Tambien se muestra que, a pesar de que la anisotropia 
es necesaria para que se generen estados entrelazados, una mayor 
anisotropia no implica necesariamente mejores condiciones para ge- 
nerar entrelazamiento que sirva para usarse en una computadora 
cuantica. 

Este trabajo esta dividido en tres partes. En la primera se ex- 
ponen nociones preliminares sobre entrelazamiento y computacion 
cuantica necesarias para comprender el estudio especifico de la ca- 
dena XY que se realiza en la segunda. La primeras dos partes estan 
divididas en dos capitulos cada una. En el Capitulo 1 se da una intro- 
duccion al concepto de entrelazamiento, una nocion del formalismo 
de la computacion cuantica y finalmente se discute la relacion entre 



VII 



VIII 



INDICE GENERAL 



estos dos. En el Capitulo 2 se presentan y discuten distintas medi- 
das de entrelazamiento. La segunda parte empieza en el Capitulo 3 
donde se presenta la cadena XY como una posible realizacion de una 
computadora cuantica, se comenta sobre las simetrias relevantes, se 
lo resuelve para dos casos simples y se muestra la capacidad de este 
sistema para generar entrelazamiento en estos casos. En el Capitulo 
4 se resuelve la evolucion temporal de la cadena XY en un caso mucho 
mas general y se realizan diversas aproximaciones que permiten des- 
cribir de manera simple su capacidad para generar entrelazamiento 
en distintos casos. Finalmente en la tercer parte se incluyen cuatro 
apendices con el objetivo de ampliar algunos temas relevantes. 

Los resultados obenidos en este trabajo ban sido aceptados para 
su publicacion en la revista Physical Review A en un trabajo conjunto 
con Raul Rossignoli. En la publicacion aprarecen ademas resultados 
que incluyen: cadenas con cantidad par e impar de qubits, resonan- 
cia de pares entrelazados de paridad positiva, un estudio detallado 
del comportamiento para v = g j b = donde la evolucion es estricta- 
mente periodica y resultados de diagonalizacion directa para n = 4. 



Parte I 

Nociones Preliminares 
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Capitulo 1 



Entrelazamiento y 
Computacion Cuantica 

La idea de hacer maquinas que calculen que se basen en las reglas 
de la Mecanica cuantica muchas veces se le atribuye a Richard Feyn- 
man. Sin embargo no fue hasta algunos anos despues de su muerte 
que, basados en los trabajos de Deutsch, independientemente Shor 
y Grover descubren algoritmos que permiten resolver mas eficiente- 
mente algunos problemas que una computadora clasica. Estos traba- 
jos fueron el puntapie inicial que llevo a armar un esquema en el que 
la Computacion Cuantica resulta una opcion viable. La computacion 
cuantica explota un aspecto fundamental de la mecanica cuantica: el 
entrelazamiento . 



1.1. Entrelazamiento y Correlaciones 
Clasicas 

El dia a dia del mundo que percibimos puede ser bien descripto 
mediante la fisica clasica. En el observamos continuamente correla- 
ciones entre eventos: el FMI impone, por medio del Gobierno local, 
una nefasta ley de educacion superior - los estudiantes manifiestan, 
los estudiantes manifiestan - sale la policia a reprimir, etc; o pongo 
un disco de Pixies en el minicomponente - escucho Where is my Mind. 
De todos estos eventos hay algunos que esta mas correlacionados que 
otros. Tomemos el ejemplo del disco de Pixies. Si tengo 3 albumes de 
esta banda pero en solo uno esta ese tema entonces hay una correla- 
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cion de 1/3 de escuchar el tema (suponiendo que el equipo anda bien y 
ninguna otra cosa sale mal). Si en cambio estoy escuchando Where is 
my Mind hay una correlacion maxima, de valor 1, de que el disco que 
puse en el minicomponente sea el de Pixies y no otro porque se que 
ese tema no lo tengo en ningun otro disco. Dos eventos clasicos estan 
maximamente correlacionados cuando la existencia de uno implica la 
del otro y viceversa. Clasicamente este limite no puede superarse. 

La mecanica cuantica es un conjunto de reglas que fueron desarro- 
lladas para describir el comportamiento de objetos microscopicos que, 
como decia Heisenberg, funciona; es decir describe adecuadamete la 
realidad. Desde sus comienzos y hasta hoy en dia la interpretacion de 
estas reglas es tema de discusion. La interpretacion de Copenhagen, 
a grande rasgos, dice: Un sistema cuantico evoluciona en una super- 
posicion de estados y cuando es medido por un objeto clasico colapsa 
a uno de los estados permitidos por el objeto clasico. Volviendo al 
caso de la musica: supongamos que yo y el minicomponente somos 
sistemas cuanticos. Esto es que podemos estar en superposiciones de 
escuchar/no-escuchar Where is my Mind y de estar/no-estar puesto el 
disco correcto. Podemos estar, si nadie nos mide, realmente en ambos 
estados al mismo tiempo y no solamente en uno u otro. Ahora, supon- 
gamos que mi novio es un objeto clasico y me pregunta -^Escuchaste 
Where is my Mind?-. Al hacerlo me esta midiendo y yo colapso a uno 
de los estados clasicos que el me permite y le digo si o no. Lo mismo 
pasa si va a revisar el minicomponente: el ve o no ve el disco, pero 
no ve una superposicion disco/no-disco porque el solo permite disco o 
no-disco pero no superposiciones. 

Siguiendo con este ejemplo y volviendo al tema de las correlacio- 
nes. Si yo y el minicomponente estuviesemos maximamente correla- 
cionados en el sentido clasico estariamos con certeza absoluta en al- 
guno de los estados "disco y escucho", "disco y no-escucho", "no-disco 
y escucho" o "no-disco y no-escucho", pero cuanticamente podemos es- 
tar con certeza absoluta en algun estado que sea una superposicion 
entre, por ejemplo, "disco y escucho" y "no-disco y no-escucho". Este 
tipo de estados esta mas correlacionado que los anteriores y no son 
admisibles en una interpretacion clasica del mundo. A este tipo de 
correlaciones que van mas alia de los llmites de la fisica clasica las 
Uamamos Entrelazamiento. 

Para que todo esto tenga sentido tiene que haber algiin objeto 
clasico que pueda medir o reconocer este tipo de superposiciones. 
Continuando con el ejemplo esto podria ser un amigo especial que no 
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pueda reconocer si estoy escuchando o no sino que reconozca si estoy 
en un estado escuchando/no-escuchando y que al ser medido colapse 
a este. El principio de incerteza aplicado a este ejemplo implica que 
no puede haber un aparato que mida simultaneamente superposicio- 
nes y estados "normales". Es decir que si mi amigo tiene como estado 
permitido la superposicion escucho/no-escucho le son prohibidos los 
estados escuchando y no escuchando por separado. Esto ultimo por 
mas raro que parezca en el contexto de este ejemplo, en fisica atomi- 
ca puede involucrar simplemente rotar un instrumento de medida. 

Esta intrigante propiedad de la mecanica cuantica fue senalada 
como una paradoja por Schrodinger, como una razon para decir que 
la mecanica cuantica no era completa por Einstein y sus "secuaces'0, 
y finalmente por muchos otros fisicos como una de las propiedades 
fundamentales de la mecanica cuantica. En 1964 Bell[ Bel871 estable- 
cio criterios matematicos estrictos para los limites que pueden tener 
las correlaciones clasicas y propuso varios experimentos que se rea- 
lizaron con exito [AGR82][ADR821, en los que se miden correlaciones 
que superan las clasicamente permitidas, demostrando que eran in- 
correctas las suposiciones en las que se basaban las observaciones 
de "los reaccionarios" Einstein, Podolsky y Rosen. El tema del gato 
de Schrodinger fue recien totalmente comprendido cuando se desa- 
rrollo la teoria de la decoherencial^. 



1.2. Computacion Cuantica 

Doy a continuacion una breve descripcion de algunos aspectos fun- 
damentales de la Computacion Cuantica que creo fundamentales pa- 
ra entender el contexto en el que fue pensado este trabajo. Para una 
exposicion mucho mas completa ver el libro Quantum Computation 
and Quantum Information de Nielsen Y Chuang flNCOOH o las Lecture 
Notes de John Preskil UPreH . 



1.2.1. Bits y qubits y compuertas 

El mundo descripto por la mecanica cuantica admite estados que 
son superposicion de los estados clasicamente permitidos. En este 



^Einsten, Podolsky, Rosen|lEPR35]| 
^Ver apendicelDl 
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sentido es natural pensar que una computadora cuantica tendra en- 
tonces como variables estados que son superposicion de los estados 
permitidos de una computadora clasica, esto es superposiciones de 
los valores que puede tener un bit: o 1. A estas nuevas variables las 
Uamamos qubits. 

Los posibles estados de un bit clasico se rotulan simplemente con 
los numeros y 1. En computacion cuantica para indicar como se en- 
cuentra un qubit se debe indicar en que estado esta el sistema cuanti- 
co que guarda esta variable. Una de las notaciones mas comunes para 
indicar estados cuanticos es la notacion de bra-kets de Dirac. En esta 
un sistema que esta en estado x se lo rotula \x). Un qubit \Q) que es 
una superposicion de los estados |0) y |1) de una computadora clasica 
se escribe como: 

|g> =a|0) + /3|l) donde \a\^ + \(5\^ = l (1.1) 

Asi como para sistemas cuanticos se generalize la nocion de bit a 
qubit, se generaliza la idea de compuertas logicas. 



Compuertas de un qubit 



La linica compuerta logica clasica de un bit no trivial que se pue- 
de definir es la NOT, que cambia el estado del bit. En computacion 
cuantica se la generaliza diciendo que actua linealmente sobre el es- 
tado. Si llamamos X a esta compuerta se tiene: 



X \Q) = X{a \{)) + (5\l)) = a\l) + (5 |0) 



(1.2) 



Si en vez de usar la notacion de bra-kets utilizamos una matricial en 
la base computacional (la de los kets |0) y |1)) la compuerta NOT se 
escribe como: 



X = 



1 

1 



de mode que X 



a 










a 



(1.3) 



Visto de esta manera es evidente que en computacion cuantica po- 
demos hacer otras compuertas de un qubit no triviales. En general 
como la compuerta, fisicamente, es la evolucion del sistema dada por 
un Hamiltoniano estas podran ser todas las matrices unitarias de 
2x2. Cualquiera de ellas que puede ser descompuesta en tres rotaci- 
nes y una fase. 
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En particular una compuerta muy relevante para la computacion 
cuantica es la Hadamard H. 



H = 



bien 



1 



V2 



1 1 

1 -1 



de modo que H 



V2 



a 


1 


' (3' 


_P _ 


" V2 


a — (3 



//(«|0) + /?|1)) = .M + ;,|0) 



72 



(1.4) 



(1.5) 



Compuertas de varies qubits 

Las compuertas de mas de un qubit clasicas suelen tomar una 
cierta cantidad de bits de entrada y dar como resultado de una ope- 
racion logica una cantidad menor de bits. Estas no son reversibles en 
el sentido que al hacer la operation hay informacion que se pierde 
y la entropla del sistema aumenta@. Las compuertas cuanticas son, 
por lo general, definidas como la extension reversible y lineal de su 
correspondiente clasico. 

En particular presento una compuerta de dos qubits relevante: la 
controlled not (CNOT) es la extension lineal de una compuerta que 
toma dos qubits y niega el segundo solo si el primero esta en el esta- 
do |1) y no cambia el primer qubit. Su representation matricial en la 
base computational es: 



UcN — 



10 
10 
1 
10 



(1.6) 



Compuertas Universales 

La relevancia de la CNOT viene de que, al igual que con la com- 
puertas clasicas AND OR y NOT puede construirse cualquier com- 
puerta clasica, con una composition finita de CNOT y compuertas de 
un qubit puede aproximarse tanto como se desee cualquier compuer- 
ta cuantica. La demostracion de lo enunciado esta didacticamente 



expuesta en el libro de Nielsen y Chuang UNCOOH donde ademas se 
senalan las referencias a los trabajos en los cuales se demostro esto 
originalmente. 



^Existen tambien logicas clasicas reversibles pero las usuales no lo son. 
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1.2.2. Algoritmos 

A groso modo hay tres cases de algoritmos eficientes en compu- 
tacion cuantica. Por un lado los basados en la transformada de Fou- 
rier discreta entre los cuales estan el famoso algoritmo de factoriza- 
cion de Shor, algoritmos que encuentran periodos de funciones, esti- 
madores de fases, etc. Por otro lado algoritmos de busqueda como el 
de Grover y finalmente algoritmos de simulacion de sistemas fisicos 
complicados. Presentare aqui la version mas simple de un algorit- 
mo del primer tipo - el algoritmo de Deutsh - que no tiene ninguna 
importancia real en el sentido de la eficiencia pero es facil de enten- 
der y muestra algunas de las ideas fundamentales0. A los tipos de 
algoritmos mencionados deben sumarseles los, para nada menos im- 
portantes, codigos de teleportacion cuantica mediante los cuales se 
pasa toda la informacion de un cierto estado a otro separado espa- 
cialmente del primero. 



Procesamiento en paralelo 



Como mencione anteriormente el primer trabajo que propuso un 
algoritmo que una computadora cuantica puede resolver mas eficien- 
temente que una clasica fue el de Deutsch. Este es en realidad un 
algoritmo algo inutil pero igualmente significante. Describo a conti- 
nuacion una version simplificada y mejorada. La idea es tratar de 
encontrar si una funcion / : {0, 1} {0, 1} es constante o no. Es decir 
si /(I) = /(O) /(I) 7^ /(O). Clasicamente se requeriria "preguntar- 
le" a / dos veces cual es su imagen. Osea calcular /(O) y luego /(I). 
Cuanticamente veremos que con solo "preguntarle" una vez basta. El 
algoritmo parte de un estado inicial de dos qubits. 

|Qo) = |0)|l) (1.7) 
Se le aplica una compuerta Hadamard a cada uno de los qubits. 



Hi ® H2 \Qo) = \Qi) = 



|0) + |1) 
V2 



|o)-|i) 

V2 



(1.8) 



Luego se aplica una compuerta U / que evalua una funci6n@ que al 
segundo qubit lo lleva a ?/ © /(x) y al primero no lo modifica (donde y 
es el valor del primer qubit y x del segundo y © es la suma modulo 2). 



*Para mas detalle y otros algoritmos ver el libro de Nielsen y Chuang P NCOOL 
^Debe entenderse que esta funcion es lineal en la superposicion de estados. 
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Uf \Qi) = IQ2) = <^ 




si /(I) = /(O) 



si /(I) ^ /(O) 



(1.9) 



Por ultimo se le aplica una compuerta Hadamard al primer qubit. 



Hi IQ2) = IQ3) = 



' ±10) 


no>-ii)] 


Si /(I) = /(O) 




no>-ii)i 

[ ^/2 J 


si /(I) ^ /(O) 



(1.10) 



Y ya esta: midiendo el primer qubit se sabe si la funcion es constan- 
te no. El argumento en contra de este procedimiento como compu- 
tadora cuantica eficiente es que, a pesar de que uno solo tiene que 
"preguntarle" a / cuanto vale solo una vez, se deben usar dos qubits 
y aplicar Hadamards antes y despues. Es cierto, este no es un circui- 
to muy eficiente pero fue el primero e ilustra claramente la idea de 
procesamiento en paralelo de una computadora cuantica. 



1.3. Entrelazamiento y Computacion 
Cuantica 

He dado hasta ahora una nocion de lo que es el entrelazamiento 
y de problemas que la computacion cuantica se propone resolver. El 
objetivo de esta seccion es formalizar la idea de entrelazamiento y 
finalmente relacionarlo con computacion cuantica. 

1.3.1. Separabilidad 

Como primera instancia en la formalizacion de las ideas sobre el 
entrelazamiento es fundamental introducir el concepto de estado se- 
parable. Un estado es no entrelazado cuando es separable y vicever- 
sa. Para hablar de entrelazamiento y separabilidad es obviamente 
necesario pensar en sistemas constituidos por mas de una parte - 
subsistemas- ya que sino no habria entre que generar entrelazamien- 
to que separar de que. Es decir algo no puede estar entrelazado con 
sigo mismo separado de si mismo. 
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Un sistema de objetos distinguiblesH es separable si su estado, 
descripto por un ket \'4j), puede ser escrito como un producto directo 
de los kets que describen cada subsistema0. De modo que un estado 
separable puede ser escrito en alguna base de la forma. 

\ip) = |mi) ^7712)2 • • • (1.11) 

= \m\m2 . • . TUn) (1.12) 

Un estado es no separable entrelazado si no puede escribirse de 
esta manera. 

Separabilidad - Matrices densidad 

La definicion anterior anterior de estados separables es valida 
solo para estados puros. Para describir estados mixtos es necesario 
introducir el formalismo de la matriz densidad. En terminos de ma- 
trices densidad un sistema de dos componentes Ay B descripto por 
una matriz densidad p es separable si [Wer89i1 

P = Y.'i-^PA®PB ; ^a>0 (1.13) 
a 

Esto implica que un estado general (puro mixto) es separable si 
puede ser escrito como una superposicion estadistica de estados pu- 
ros separables. Es decir si existe alguna base en la cual puede ser 
escrito de esta manera. 

1.3.2. Estados de Bell 

Un caso muy importante de estados entrelazados son los que se 
llaman estados de Bell, que a veces se esconden bajo los seudonimos 
de pares estados EPRH. Los estados de Bell son estados de un siste- 
ma compuesto por dos subsistemas -sistema bipartito- de dimension 

^Siempre consideraremos que los objetos que se describen son distinguibles. Esta suposicion 
tiene sentido en el marco de la computacion cuantica donde se requiere poder distinguir un qubit de 
otro. Por mas que usualmente se piensa en usar algiin grade de libertad de particulas identicas como 
qubit se considera que cada una esta confinada espacialmente. Las definicion de estado separable 
para un sistema de particulas indistinguibles es distinta a la que se da aqul. 

^El subsistema i en el estado es descripto por el ket |mj) .. En el caso de qubits los pueden 
valer solo o 1. 

^Nombre debido a los "reaccionarios" Einstein, Podolsky y Rosen 
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dos. Estos son: 

. Mijli) (1.14a) 

\M - (1.14b) 

|Ao> = M_Jli> (1.140 

\M = (i,i4d) 

Los estados de Bell son de gran importancia ya que, como veremos en 
la secci6n l2?l|. estos son maximamente entrelazados en el sentido que 



en un sistema bipartito de dimension dos no hay estados que tengan 
mayor entrelazamiento que estos. Ademas son cuatro y son ortonor- 
males de modo que son una buena base para este sistema. Por ultimo 
menciono que se los usa como unidad de medida de entrelazamiento: 
un estado de Bell tiene un ebit de entrelazamiento. 



Estados de Bell y compuertas 

Es facil generar estados entrelazados de Bell a partir estados no 
entrelazados usando las compuertas descriptas en la seccion 11.21 Por 



ejemplo si partimos del estado |11) y aplicamos una Hadamard al 
primer qubit se obtiene (|0) — |1)) |0) que todavia es separable. Si 
a continuacion se aplica una CNOT se obtiene (|01) — \10))\/2 que es 
el estado entrelazado de Bell \/3ii). De este mismo modo, cambiando 
los estados iniciales, se pueden generar todos los estados de Bell. 

Del procedimiento anterior se ve que quien realmente "genera" el 
entrelazamiento es la compuerta CNOT. Sin embargo la aplicacion 
de la compuerta CNOT directamente sobre el estado lo lleva al 
1 10) que tampoco es entrelazado. Volviendo a las ideas de la intro- 
duccion esto puede explicarse de la siguiente manera. El estado |11) 
tiene correlaciones clasicas entre el valor de cada qubit. La compuer- 
ta CNOT sobre este cambia las correlaciones entre un qubit y el otro. 
Es decir, al alimentar la CNOT con estados de la base computacional 
esta los hace interactuar como sistemas clasicos y no genera entrela- 
zamiento. En cambio si alimentamos a la CNOT con un estado que es 
superposicion de estados en la base computacional esta reconoce las 
caracteristicas cuanticas del sistema y, al hacerlos interactuar, gene- 
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ra entrelazamiento. El entrelazamiento se genera en la interaccion 
de sistemas cuanticos. 



1.3.3. Entrelazamiento y Computacion Cuantica 

El algoritmo de procesamiento paralelo presentado en la seccion 
L2] puede parecer un ejemplo bastante malo a los fines de mostrar 



como la computacion cuantica usa el entrelazamiento como recurso. 
Explicare ahora porque no es tan asl y de este modo trato de dar una 
idea del rol que cumple el entelazamiento en los algoritmos de compu- 
tacion cuantica. Las bases de la observacion mencionada se basa en 
que en todos los pasos del algoritmo los estados son separables (ver 
ecuaciones II. 7|. II. 8|. II. 9|. II.IOD . Sin embargo examinando la forma que 



tiene Uf por ejemplo para el caso en que /(I) = 1 y /(O) = se en- 
cuentra que es exactamente una compuerta CNOT y que en los otros 
casos no es mas que composiciones CNOTs con Hadamards. Mas aun 
los estados a los que se les aplica [// no le son propios de modo que 
Uf los hace interactuar como objetos cuanticos y no clasicos. A pesar 
de que no se genera entrelazamiento en ninguna etapa del proceso la 
interaccion entre los qubits es de origen cuantico y es esta la que per- 
mite realizar la tarea del modo que se hace. Es mas, es posible que 
este algoritmo no sea realmente eficiente debido a no hay estados 
entrelazados de por medio. 

Finalmente hago enfasis en que cualquier algoritmo puede des- 
componerse en CNOTs y compuertas de un qubit y, en particular, 
practicamente todos los algoritmos eficientes de computacion cuanti- 
ca usan los CNOTs para generar estados entrelazados. El rol que 
cumple el entrelazamiento en la computacion cuantica ha sido muy 
discutido en particular en torno a las realizaciones que se basan en 
tecnicas de NMR -resonancia magnetica nuclear- [iLPOlJ . Contraria- 
mente el rol del entrelazamiento en los codigos de teleportacion es 
fundamental. 



Capitulo 2 

Medidas de 
entrelazamiento 



^Que se espera de una medida de entrelazamiento? 

Una medida de correlaciones puramente cuanticas. Es decir que dis- 
tinga ente correlaciones cuanticas y clasicas y solo considere las de 
origen cuantico. Una medida que distinga entre estados separables 
y no separables. En este camino se ban dado grandes avances en los 
ultimos 10 anos. No se ba llegado a encontrar una medida universal 
de entrelazamiento aun. Por ejemplo Peres[Per96] encontro un crite- 
rio de separabilidad que es necesario pero no suficiente. Se ban pro- 
puesto varias medidas aditivas y no aditivas DGGOU . En este trabajo 
me concentrare en dos medidas aditivas: la entropia de non Neumann 
que sirve solo para estados puros, y el entrelazamiento de formacion 
que sirve para estados mixtos [IBDSW96ll y coincide con la anterior 
para estados puros. A pesar de que el entrelazamiento de formacion 
es muy dificil de calcular para un caso general Hill y Wooters [ HW97II 
encontraron una forma explicita de calcularla para un par de qubits 
como funcion de una cantidad que llamaron Concurrencia. En pos 
de unificar medidas simples de entrelazamiento damos aqm un paso 
mas proponiendo una medida equivalente a la entropia de un qubit 
que sea comparable a la Concurrencia. 



2.1. Entropia - Estados Puros 

La entropia de von Neumann de un subsistema es una buena me- 
dida de entrelazamiento para estados puros. Mide el entrelazamiento 
entre dos subsistemas sin importar su dimension. Daremos a conti- 



13 



14 



Capitulo 2. Medidas de entrelazamiento 



nuacion la definicion, luego justificaremos porque esta es adecuada y 
finalmente estudiaremos algunos ejemplos. 

2.1.1. Definiciones 

La entropia de von Neumann adimensional0 de un sistema cuanti- 
co caracterizado por una matriz densidad p se define como: 

S{p) = -Tr[p\og^p\ (2.1) 

La matriz densidad reducida de un subsistema se define como traza 
parcial sobre todas las variables que no pertenecen al subsistema. 

Psuh = Tr^.^iXp] (2.2) 

Supongamos un sistema bipartito, que puede ser dividido en dos sub- 
sistemas disjuntos A y B. Toda la informacion del estado de cada 
subsistema esta en las matrices densidad reducidas {pA y p^) de cada 
subsistema. 

PA = TtbIpab] ; Pb = TtaIpab] (2.3) 

La entropia de un subsistema A se define analogamente a la entropia 
para el sistema entero pero con la matriz reducida pA en vez de la 
matriz densidad entera p. 

S {p a) = -Tr[p Alog2 pa] ; S{pb) = -Tr[pBlog2pB] (2.4) 

Ademas es inmediato ver, usando las propiedades de la descomposi- 
cion de Schmidt (el apendice 13) que los autovalores no nulos de pA y 
Pb son iguales de modo que sus entropias son iguales S{pa) = S{pb). 
Se puede definir una buena medida de entrelazamiento E{pab) entre 
dos subsistemas Ay B, que forman un estado puro, como la entropia 
de cualquiera de los subsistemas. 

E{pab) = SipA) = S{pb) (2.5) 



^En el contexto de Teoria de la Informacion Cuantica se la suele definir adimensional y en base 
2 para que, como medida de desorden, valga 1 para un qubit completamente desordenado y para 
uno completamente ordenado. 
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2.1.2. Propiedades/Justificacion 

Doy a continuacion una serie de propiedades de la entropia de 
un subsistema, que justifican su eleccion como medida de entrelaza- 
miento entre subsistemas de un estado puro. Algunas de estas las 



ideas salen de un resumen hecho por Bennett et. al. [IBDSW96,] 



El entrelazamiento de A con B es igual al de B con A . Esto se 
desprende directamente de la definicion que considera un sistema 
bipartite puro y de la observacion hecha anteriormente sobre los 
autovalores de los subsistemas. Osea S{pa) = S{pb] 



Para estados separables E{pab) = . Si un estado puro es separa- 
ble pab = Pa® Pb ,donde paY Pb son estados puros, de modo que 
E{pab) = S{pa) = S{pb) = dado que la entropia de un estado 
puro es cero (ver apendice |Al). 

Los estados de Bell estan maximamente entrelazados. Un cal- 

culo directo conduce a E{\/3ij) {/3ij\) = 1, Vi, j 

Es aditiva E{pi ® p2® ■ ■■pn) = E{pi) + E{p2) + . . . E{pn). El entrela- 
zamiento de sistemas un conjunto de sistemas independientes es 
la suma del entrelazamiento de cada sistemaE 

E{pab) se conserva frente a operaciones locales unitarias . 

Transformaciones locales del tipo U = Ua ® Ub no cambian el 
entrelazamiento entre los subsistemas [IBBPS96II . 



E{pab) no se incrementa con operaciones locales no unitarias. 

Transformaciones locales no unitarias, como mediciones, solo pue- 
den disminuir el entrelazamiento pero no incrementarlo. 

El entrelazamiento puede diluirse o concentrarse@ Esto es: 
realizando solamente oparaciones locales en cada subsistema y 
permitiendo comunicacion clasica entre ellos se puede, 

a partir de n copias de pab, preparar m sistemas identicos p'^^ 
con eficiencia m/n « E{pAB)/E{p'ji^^) en el sentido asintotico de n 
grande.[BBPS96Jl 



^Algunos autores llGGOllI argumentan que no es necesario exigir aditividad para una medida de 
entrelazamiento. 

^Esta propiedad es la base teorica utilizada para construir otras medidas de entrelazamiento 
como entropia de formacion y la herramienta fundamental de diversos codigos correccion de errores. 
*LOCC - Por sus siglas en ingles Local Operations and Classical Communication 
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2.1.3. Entropia de un qubit 

La entropia de un qubit es una medida del entrelazamiento de ese 
qubit con el resto del sistema, por eso la llamare entrelazamiento de 
uno con el resto. 

La matriz reducida de un qubit siempre puede escribirse en la base 
computacional como: 

a b 
6* 1-a 



Pi 



(2.6) 



y tiene por autovalores A± = | (l =t — 2ay + |26|2^ de modo que 
la entropia de un subsistema arbitrario de un qubit es 

S{pi) = -A+logA+ - A_logA_ 

= -A+logA+-(l-A+)log(l-A+) (2.7) 

En el caso que el sistema completo tenga paridad de espm definida 
(ver apendice [B]) se tiene que 6 = y la entropia solo depende solo de 
a y es simplemente 

= -a log a -{I -a) log(l - a) (2.8) 

Ademas a = (Sz) + 1/2 de modo que si pi conmuta con P alcanza con 
conocer {Sz) para obtener una medida de la entropia. 



2.2. Entrelazamiento de Formacion - Es- 
tados Mixtos 

El entrelazamiento de Formacion es una buena medida de entre- 



lazamiento para estados mixtos [iBDSW96ll . Mide el entrelazamiento 



entre dos subsistemas sin importar su dimension. Tiene la inconve- 
niencia de ser muy diflcil de calcular. Es una extension de la entropia, 
como medida de entrelazamiento, para estados mixtos. 
La idea detras de la definicion se basa en los mecanismos de desti- 
lacion y dilucion mencionados en la seccion anterior y la idea egip- 
cia elegir una unidad estandar de referenda, as! como el kilogramo 
el metro patron. Se toma como unidad estandar de referenda de 
entrelazamiento el entrelazamiento de un singlete puro. Por meca- 
nismos de destilacion se pueden generar n singletes puros a partir de 
m copias de un estado mixto o bien diluyendo se generan m estados 
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mixtos a partir de n singletes puros. El entrelazamiento de formacion 
de un estado es el cociente m/n en el sentido asintotico de n grande. 



2.2.1. Definicion 

Para un sistema bipartito puro el entrelazamiento de formacion 
E{pab) es la entropia de von Neumann de uno de los subsistemas. 

E{pab) = S{pa) = S{pb) (2.9) 

Para un ensamble de sistemas bipartitos puros S = {pi,pAB} 
trelazamiento de formacion E{S) se define como el promedio de los 
entrelazamientos de formacion de cada estado /o^^ pesadas con la 
probabilidad pi del estado. 

E{S) = J2pMp'ab) (2.10) 

i 

Para un sistema bipartito mixto pab el entrelazamiento de forma- 
cion E{pab) se define como el mfnimo de los E{S) sobre todos los 
posibles ensambles S = {pi,pAB} representen al estado mixto 

PAB = J2^P^PAB^ 

E{pAB)=^i^Y.P^^^PAB) (2.11) 



2.2.2. Propiedades/Justificacion 

En el trabajo de Bennett et. al. []BDSW96 1 en el que proponen al 
entrelazamiento de formacion como medida de entrelazamiento de- 
muestran que para generar un estado bipartito pab solo a partir de 
operaciones locales y comunicacion clasica entre los subsistemas es- 
tos deben compartir el equivalente a E{pab) singletes puros. El en- 
trelazamiento de formacion E{pab) es la cantidad de entrelazamien- 
to necesaria para formar pab- Ademas demuestran que E{pab) no 
puede incrementarse con ningun tipo de operaciones locales y comu- 
nicacion clasica. Es decir el entrelazamiento de formacion tiene las 



mismas propiedades (ver seccion 12. 1.21) que la entropia pero para es- 
tados mixtos. 
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2.3. Concurrencia 

La concurrencia es una medida de entrelazamiento derivada a 
partir del Entrelazamiento de Formacion entre un par de subsiste- 
mas de dimension 2 cada uno. Tiene la bondad de que es facil de 
calcular. Fue originalmente propuesta por S. Hill y W.K. Wootters 
^HW97il y refinada por W.K. Wooters nWoo98ll . 



2.3.1. Definicion 

La concurrencia C de un sistema arbitrario de dos qubits i,j se 
define como: 

C = Max[2Xm - TrR, 0] ; R = \pTppT\^''^ (2-12) 

Donde = Tr^{ijj[p] (traza parcial sobre todas las variables que no 
son z, j); Xm es la maximo autovalor de R;yp2 = 4:S^SjP2SfSj (matriz 
densidad p2 con los espines invertidos). 

A los efectos del calculo de la concurrencia es equivalente usar la 
raiz de los autovalores de la matriz no hermitica pp2 ya que las cuen- 
tas son mas simples y los resultados los mismos. En este sentido, si 
Ai, A2, A3, A4 son las raices de estos autovalores en orden decreciente 
la concurrencia es simplemente@ 

C = Max[Xi -X2-X3- A4, 0] (2.13) 



2.3.2. Propiedades/Justificacion 

En los trabajos P HW97ll . nWoo98ll Hill y Wootters demuestran que 
la entropia de formacion para un estado arbitrario de 2 qubits es 
una funcion monotonamente creciente de la concurrencia y por ende 
proponen usarla como medida de entrelazamiento. Explicitamente el 
resultado que obtienen es que la entropia de formacion es 

E(M) = H (i±4^) (2.14) 

h{a) = — a log a — (1 — a) log(l — a) 



^Esta ecuacion tambien vale para los autovalores de R. 
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2.3.3. Cadena ciclica con paridad de espm 

Mas adelante estudiaremos una cadena ciclica cuyo estado tiene 
paridad definida. En este caso la matriz densidad de dos qubits en la 
base computacional se escribe@: 



P2 = 



/ a d*\ 

6 e 

e 6 

\d c J 



(2.15) 



Un calculo directo conduce al resultado: 



C = 2 Max[\e\ - ^Jac, \d\ - 6, 0] 



(2.16) 



Donde solo uno de los valores puede ser positivo ya que la positividad 
de p2 implica que |e| < b y \d\ < ^/ac. 

Es importante tambien tener en cuenta como reconstruir a partir 
de los valores medios de los operadores relevantes. En este sentido se 
calculan directamente como (O) = Tr\Op'^ los valores medios. 



ha-2b + c) 



= d 



(2.17) 



La ecuacion que falta sale de la condicion Tr[p2\ = 1; de este modo 
1 = a + 26 + c. Estas son todas las herramientas necesarias para 
reconstruir p2 y de este modo la concurrencia de pares. 



2.3.4. Concurrencia de estados de Bell 

Los estados de Bell tienen todos maxima concurrencia C = 1. Para 
los estados mixtos que conmutan con la paridad de espfn(ver apendi- 



ceEl) vemos que si el termino de 12. 161 que da a lugar a la concurrencia 



es — 6 entonces el entrelazamiento corresponde a un estado par 
mientras que si el termino relevante es |e| — ^/ac corresponde a uno 
impar. La paridad de espm resulta muy util para clasificar los esta- 
dos y el tipo de entrelazamiento de pares en los sistemas a estudiar. 



^Que haya terminos iguales en la diagonal de la matriz es debido a la simetria ciclica de la 
cadena, los terminos nulos son debido a la simetria de paridad (ver apendicelBll. 
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0.2 0.4 0.6 0.8 

a(t) 



Figura 2.1: Concurrencia 2^ a{\ — a) (Imea solida) y entrelazamiento de forma- 
cion h{a) = —a log a — (1 — a) log(l — a) (linea punteada) como funcion de a 



2.3.5. Concurrencia de uno con el Resto 

Tanto el entrelazamiento de formacion de un par como el de uno 
con el resto se escriben como: 

h{a) = -alog a- {I- a) log{l- a) (2.18) 

donde para el entrelazamiento de un par a = (1 + Vl — C'^)/2 y para 
el entrelazamiento de uno con el resto a = A+ o a = A . El entrela- 
zamiento de formacion de un par se reduce a la entropia de un qubit 
en el caso en que el par es un estado puro. En este sentido ambas 
medidas son comparables pero tienen la desventaja de que haya lo- 
garitmos de por medio. 

Debido a que C = 2 ^/a{l — a) y que como fue senalado anteriormen- 
te C y h{a) son funciones monotonamente crecientes una de la otra 
puede definirse la concurrencia de uno con el resto de este modo pero 
usando a = A+. 



Ci(M) = 2^ya{l-a) = 2^/X+{l- X+) (2.19) 
C2(M) = 2^a{l-a) = C (2.20) 

Donde hemos acunado los sfmbolos Ci y C2 para referirnos a las 
medidas de entrelazamiento de uno con el resto y de pares respec- 
tivamente. En lo que continua usare los terminos entrelazamiento 
y concurrencia indiscriminadamente para referirme siempre a estas 
medidas. 

Un ultimo comentario que merece la atencion es la forma de la 



funcion 2y^a(l — a). En la figura [2J] se muestra esta funcion como la 
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alternativa ah{a). Ambas estan definidas en el intervalo < a < 1 y 
tienen un maximo en a = 1/2 donde valen 1. 
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Parte II 

Entrelazamiento en una 

Cadena XY 
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Computacion cuantica 
en una cadena XY 



De los muchos posibles sistemas que se podrian elegir para ha- 
cer computacion cuantica, y sobre los cuales se esta investigando ac- 
tualmente [|roa02][NC 001 , algunos pueden ser descriptos aproxima- 



damente con Hamiltonianos XY. En particular estos podrian ser de 
NMR(resonancia magnletica nuclear), de Estado Solido (basados en 
espfn como en "Quantum dots") o superconductores. En este sentido 
se ha estudiado el entrelazamiento de estados termicos para este ti- 



po de cadenas URCOSII . Tambien se ban propuesto esquemas mediante 
los cuales se podrian implementar compuertas cuanticas en cadenas 
de este tipo de sistemas[BB03]. 

3.0.6. El Hamiltoniano 

El sistema que estudie es una cadena clclica de n espines inter- 
actuando por medio de acoplamientos de Heisenberg del tipo XY a 
primeros vecinos en presencia de un campo magnetico uniforme y 
constante perpendicular a la direccion de las interacciones. El Ha- 
miltoniano para este sistema es 

n 
i=l 

donde n + 1 = 1. Definiendo los operadores usuales de subida 5+ y de 
bajada 5 como: 

s+ = s^ + isy ; s; = s^^-isy 0.2) 
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el Hamiltoniano queda expresado como: 

n 

n 

= E K - i^SrS-,, + gSrSt, + h.c.)] (3.3) 

i=l 

Donde {v, g) = {v^ ± Vy) /4. Kvlo llamare parametro de "hopping" y a 
g anisotropia. 



3.0.7. La condicion inicial 

En todos los casos estudiare la evolucion temporal del entrelaza- 
miento de un sistema que inicialmente tiene todos los espines alinea- 
dos antiparalelos al campo magnetico. Tiene sentido estudiar esta 
condicion inicial ya que es experimentalmente muy facil de generar. 
Un campo magnetico suficientemente fuerte siempre va a alinear los 
espines en una misma direccion. La funcion de onda a t = es 

|*(0)) = Ui...i) (3.4) 
el estado del sistema como funcion de t es (se asume en lo sucesivo 

h = i) 

|^(t)) = e-^^^Ui...i) (3.5) 



3.0.8. Paridad de espm 
Paridad y Hamiltoniano 

Es simple ver que Hamiltoniano conmuta con el operador paridad 



definido en el apendice IB. 2 



Recordemos que este operador distingue la paridad de un estado si 
hay una cantidad par o impar de espines para "arriba". El Hamil- 
toniano tiene terminos es evidente que conmutan con P. Tiene 
tambien terminos tipo S^S^_^^ y S^S^_^^. Los primeros tampoco cam- 
bian la paridad ya que "suben" un espm y "bajan" otro sin alterar la 
cantidad de espines que hay "arriba". Osea, sin modificar la paridad. 



3.1. Cadenas XY Cortas 
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Los segundos "suben" dos espines por lo que tampoco cambia la pa- 
ridad. Como ningun termino del Hamiltoniano cambia la paridad de 
una funcion de onda dada se tiene que: 

[P, H] = (3.7) 



Paridad y la evolucion temporal 

El estado inicial conmuta con P (tiene paridad definida). Esto se 
ve directamente de su definicion 13.41 Como P conmuta con H y con 
p(0) se obtiene inmediatamente que p{t) tiene paridad definida para 
cualquier t. 



3.1. Diagonalizacion directa: 

Entrelazamiento en cadenas cortas 



3.1.1. Dos qubits 

Empezaremos estudiando el caso simple de 2 qubits. Por las razo- 
nes mencionadas en la seccion 13.0.81 el sistema evolucionara dentro 
del subespacio de paridad definida que contiene al estado inicial \ 
Este es el subespacio de paridad positiva. En la base normal de auto- 
vectores de (|tt) y Ijj)) el Hamiltoniano en forma matricial es: 



H 




(3.8) 



que tiene autovalores E± = ±X con 



(3.9) 



de modo que la evolucion temporal del sistema, dada por la ecuacion 
3M es 



|*(t)) = z|sin(At)|TT) + 



cos(Ai) + i— sm{Xt) 
A 



III) 



(3.10) 



En este caso al haber solo dos qubits el entrelazamiento de pares es el 
mismo que el de uno con el resto: Ci = C2. Para calcularlo es necesario 
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conocer la funcion a{t) (ver seccion 12. 3. 51) . Esta es simplemente una 
lorenziana modulada con centro en 6 = 0. 



a{t) 



y su maximo se da para t 



9 



52 + (/^ 

(2m+l)7r 
2A 



, m eTZ 

9' 



(3.11) 



52 + ^2 



(3.12) 



El estado maximamente entrelazado de dos qubits (se da cuando 
a{t) = 1/2. Si am > 1/2 siempre se logra un estado de Bell para 
tm = arcsin((2aTO)~^/2); este caso se da cuando \b\/g < 1. Si < 1/2 
nunca se llega a generar un estado de Bell perfecto (sin ruido), el en- 
trelazamiento maximo se da para tm = 7r/2A cuando a{t) = a^. La 

alcanza- 



figura [371 . II muestra el entrelazamiento maximo 2y^a(l — a) 
do en funcion de b/g. Se observa una meseta estricta -plateau- para 
\b\/g < 1 y un comportamiento asintotico para \b\/g 1 de la. forma 
2g/\b\. El entelazamiento en este caso no depende del parametro de 
"hopping" V sino que solo depende de la relacion entre la anisotropia 
de la interaccion g y e\ campo magnetico b. 




Figura 3.1: Entrelazamiento Ci,2 maximo como funcion de b/g para dos qubits. 



Veamos las caracteristicas mas importantes de la evolucion tem- 
poral para distintos rangos de b/g. De la ecuacion I3.10|. que rige la 
evolucion temporal de este sistema y de la condicion a{t) = 1/2 se ve 
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que los estados maximamente entrelazados que se pueden obtener 
son 



V2 



(|TT) + e^^1ii)) 



cos <p = - 
9 



(3.13) 



Este estado no es ninguno de los estados de Bell en la base normal 



(11.141) pero tiene el mismo entrelazamiento que un estado normal de 
Bell y tiene paridad positiva, no es mas que un estado de Bell rotado 
y con una fase. Ademas se ve claramente, y puede verificarse con las 



figuras 13. 2L que para b/g :$> 1 el termino relevante es el que tiene al 
estado inicial |||) y por lo tanto el entrelazamiento es pequeno; mien- 
tras que para b/g < 1 los terminos relevantes son los que tienen los 
senos, que multiplican tanto a Ijj) como a ItT). de modo que efectiva- 
mente se genera un estado de Bell. 




(a) b = 0,2 



(c) 6 = 2 





(b) 6 = 1 




(d) 6 = 4 



Figura 3.2: Entrelazamiento Ci 2 en funcion del tiempo en unidades de h/v para 
dos qubits y distintos valores de b/g. 



Una particularidad interesante para notar, en vistas de aplicar 
esto a sistemas reales que funcionen como computadoras cuanticas, 
es que el tiempo que permanece el sistema con Ci 2 ~ 1 es mucho 
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mayor cuando h/g = 1 (fig |3.2(b)D . Es decir que no por agrandar la 
anisotropia se tienen mejores condiciones de formacion de entrelaza- 
miento ya que para h/g < 1 el tiempo que permanece el sistema en 
Ci 2 ~ 1 es menor (fig |3.2(a)D . En sistemas de "Quantum dots", por 
ejemplo, se ban reportado llIAB+9 9] valores de g ^ 0,lmeV que impli- 
can una escala de tiempo del orden de las decenas de pico segundos 
(10-lls). 



3.1.2. Tres qubits 

Seguimos estudiando el caso un poco mas complejo de 3 qubits. En 
este caso el estado inicial es por lo que el subespacio relevante 
es el generado por los vectores |Ui) ,|iTT),|TTi) y iTiT)- De los ultimos 
tres estados, debido a la simetria ciclica del Hamiltoniano, solo se 
puebla la combinacion simetrica de ellos y por lo tanto el espacio re- 
levante tiene, todavia, dimension 2. Elegimos para este subespacio la 
base que contiene al estado inicial y al estado de Werner [Wer89l 
IW) = (UTT) + iTTi) + ITiT)) /V3. El Hamiltoniano es en esta base 

" - ( ^-t ) 

que tiene autovalores E± = Eq±\ con 

E^ = -^--v ■ X=y/{b-vy^3g^ (3.15) 

de modo que la evolucion temporal del sistema, que viene dada por la 
ecuaci6n l3.5|. es 



/ X N .b — V . , - - 
cos[Xt) + I — - — sm[Xt) 
A 



|TTT) + ^sin(A^)|W) (3.16) 
A 



En este caso, al haber tres qubits, hay diferencia entre el entrelaza- 
miento de uno con el resto y de pares, sin embargo veremos que estan 
mtimamente relacionados y que los resultados que se obtienen pue- 
den entenderse muy bien en terminos intuitivos de entrelazamiento. 
La funcion a{t) es, como antes, una lorenziana modulada pero con 
centro enb = v, ancho g y valor maximo 2/3. 

«W = n — ^^2' o 2 sin'(^^) (3.17) 

[b — vy + Sg"^ 
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y su maximo se da para t = m e TZ 

«•» = (6 - v)2 + 3g2 <3.18) 

Haciendo un analisis similar al que se hizo para dos qubits (seccion 
13.1. ID se ve que > 1/2 para (6 — / < 1 y por ende existe un 
tiempo itm = A"^ arcsin((2a„^)~^/^)) para el cual se llega a un estado 
maximamente entrelazado. En cambio cuando {h—v)'^/g^ > 1 nunca se 
a llega generar un estado maximamente entrelazado dado que a{t) < 
1/2 para todo t. El comportamiento del entrelazamiento de uno con el 
resto es por lo tanto igual que en el caso de dos qubits (fig 13.1. ID pero 
la meseta esta desplazada en v (limites dados por (6 — /g^ < 1) y 
cae asintoticamente para h > {v, g) como Ci « 2^/2g/h. 




Figura 3.3: Entrelazamiento de pares C2 como funcion de a{t) para tres qubits. 



El entrelazamiento de pares, al tratarse de mas de dos qubits, 
difiere del de uno con el resto. Expresado como funcion de a{t) el en- 
trelazamiento de pares para este sistema es 



Co = 



a{t) - V[2 - 'ia{t)]a{t) 



(3.19) 



La figura 13^ muestra esta funcion para < a{t) < 2/3, que es el 
rango de valores que puede tomar a{t). Observese primero que nunca 
llegan a haber estados maximamente entrelazados C2 < 1 Va(t). Hay 
dos regiones que se distinguen por la anulacion de C2 en a{t) = 1/2. 
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Para a{t) por debajo de este Ifmite los estados de Bell responsables 



del entrelazamiento (ver seccion 12.31) son los que tienen la misma pa- 
ridad que el estado inicial y el maximo de esta parte se da cuando 
a{t) = 1/6 donde C2 = 1/3. Para valores de a{t) > 1/2 hay tiempos 
en los cuales se generan estados entrelazados con paridad opuesta 
a la del estado inicial, aquf el maximo se da para a{t) = 2/3 donde 
C2 = 2/3 que es justamente el entrelazamiento de pares del estado 
de Werner. Es posible generar estos estados de paridad opuesta solo 
cuando se esta cerca de la resonancia, es decir cuando (6 — /g'^ < 1. 



En la figura 13.41 se muestran los entrelazamientos mencionados. El 
entrelazamiento de uno con el resto, que tiene la misma forma que 
para el caso de dos qubits, el entrelazamiento de pares con paridad 
positiva que tiene un una meseta para \{b — v)/g\ < 3 y luego cae 
asintoticamente para b > {v,g) como C2 ~ \/2g/\b\ y el entrelaza- 
miento de pares con paridad negativa que tiene un pico pronunciado 
centrado enb = v de altura 2/3 y se anula para |(6 — v)/g\ > 1. 




Figura 3.4: Entrelazamiento maximo de uno con el resto Ci (Imea solida), de 
pares C2 con paridad positiva (linea a trazos largos) y de pares C2 con paridad 
negativa (linea a trazos cortos) como funcion deb/v para g = I. 



A continuacion mostramos la evolucion temporal (fig l3.5D para ca- 
da una de las zonas. En la resonancia (fig |3.5(a)D se ven las carac- 
teristicas mas relevantes del sistema: los picos de la concurrencia que 
corresponden a estados con paridad positiva (picos suaves) y negativa 
(picos agudos) y el entrelazamiento de uno con el resto que llega a la 
saturacion en tiempos distintos a los que aparecen los picos de la con- 
currencia. En el punto critico (fig |3.5(b)| ) , es decir donde (6 — v)/g= 1, 
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el entrelazamiento "se queda" saturado por un tiempo (como en el 
caso de dos qubits) pero desaparece el pico de la concurrencia para 
estados con paridad negativa. A medida que los parametros se alejan 



de la resonancia los picos de la concurrencia colapsan (fig |3.5(c)D y 
finalmente el comportamiento es similar al caso de dos qubits, el en- 



trelazamiento de uno con el resto sigue a la concurrencia (fig |3.5(d)D , 
pero siempre la concurrencia es menor. 





(a) (b - v)/g = 



(b) {b - v)/g - 1 





(c) (6 - v)/g = 2 



(d) (b - v)/g = 5 



Figura 3.5: Entrelazamiento de uno con el resto (lineas solidas) y entre pares 
(Imeas punteadas) en funcion del tiempo en unidades de h/v para tres qubits y 
distintos valores de {b ~v)/g. 
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Capitulo 4 

Transformaciones 
Canonicas: 
Entrelazamiento en 
cadenas XY largas 



4.0.3. Introduccion 

En este capitulo resolvemos el problema de la evolucion temporal 
de la concurrencia de pares y de uno con el resto para el sistema 
descripto por el Hamiltoniano del capitulo [3] para una cantidad impar 
arbitraria de qubits. 



El procedimientor LSM61ll consiste en hacer una transformada de 
Jordan- Wigner [iJW28(l que mapea los operadores de espin S^, a 
operadores puramente fermionicos c-,q y luego diagonalizar el Ha- 
miltoniano mediante una Transformacion Canonica, que conserva 
las relaciones de conmutacion. Una transformacion canonica puede 



ser descompuesta en tres partes nRS41ll : una transformacion unitaria 



que mapea los operadores q en operadores di del tipo di = Ylij MijCi, 
una trasformada especial de Bogoliubov (o transformada BCS) que 
mapea los operadores di en operadores de cuasiparticulas y una 
ultima transformacion unitaria entre los operadores a^. La primer 
transformacion unitaria la descompongo en dos partes: un cambio de 
fase dependiente de la paridad elegida para eliminar el termino no 
cuadratico que introduce la transformada de Jordan- Wigner con con- 
diciones ciclicas y una transformada de Fourier discreta que lleva al 
Hamiltoniano a una forma canonica. Finalmente hago la transfor- 
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mada BCS que deja el Hamiltoniano en una forma diagonal de modo 
que la ultima transformacion unitaria no es necesaria. Una vez dia- 
gonalizado el Hamiltoniano calculo exactamente los valores medios 
en funcion del tiempo necesarios para dar las formas explicitas de 
las concurrencias (ver ecuaciones I2.17D . 

Para el caso de tres qubits los resultados obtenidos coinciden exacta- 
mente con los de la secci6n [3.1.2l 

Finalmente se hace un analisis de la dependencia funcional del maxi- 
mo de las concurrencias con el campo magnetico b para distintos ran- 
gos de v,g (parametro de "hopping", anisotropia de la interaccion). 



4. 1 . Diagonalizacion 

4.1.1. Transformacion de Jordan- Wigner 

Partimos del mismo Hamiltoniano de espines 1/2 en un campo 
magnetico que usamos en la seccion anterior. 

n 

^ = E [^^^ - + gsrst, + h.c.)] (4.1) 

Las relaciones de conmutacion de los oper adores de subida y bajada 
(S:^ , 5^ ) son en parte fermionicas 

{5-,5+} = l; {S;f = {S-f = (4.2) 

y en parte bosonicas 

[S^, SJ] = K+, 3+] = [Sr, Sr]=0;zj^j (4.3) 

Para poder diagonalizar el Hamiltoniano con transformaciones canoni- 
cas (Fourier y Bogoliubov) como haremos luego es preciso que las re- 
laciones de conmutacion de los operadores sean fermionicos o bosoni- 
cas, pero no una mezcla como es el caso de , 5^ . Para llevar estos 
operadores a una forma puramente fermionica hacemos ahora una 
transformacion de Jordan- Wigner [JW281 que mapea los operadores 
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y S en operadores cJ,Cj puramente fermionicos. 

i-i i-i 
Cj = exp [iri S^S^]S[' Slj^ = exp[—7ri''^^c^jCj]ci (4.4) 
i=i 3=1 

i-l i-1 

c] = S^exp[—TTi ^ 'S'l'^ = c]exp[7Ti ^ cjc-,] (4.5) 

donde se tiene que: 

{q, cj} = (5^^-, {c„ Cj} = {c|, cj} = (4.6) 

Para poder escribir el Hamiltoniano en terminos de estos nuevos 
operadores vemos que de las definiciones de las transformaciones se 
desprende que: 

cjc, = S+Sr (4.7) 

Ademas, dado que cjc^ es un numero de ocupacion que puede ser o 
1 se tiene que: 

exp[7riclci] = exp[—7riclci] (4.8) 
Ademas para i = 1,2, . . . , N — 1 se tiene que: 

cjci+i = Sl^S-^i ; c-c-_^-^ = SlSl_^i (4.9) 

de modo que, para el caso de extremos libres, el Hamiltoniano es 

n—l ^ 

^ = E - 2) - + + h.c.)] (4.10) 

El caso que nos interesa resolver, sin embargo, es el de la cadena 
ciclica. Debemos entonces tener en cuenta los terminos 

S+S^ =-4ciP T^cl^ci (4.11a) 
S+St =-cic\P ^cic\ (4.11b) 

donde P es el operador paridad de espfn (apendice E]). De este modo 
se obtiene el Hamiltoniano para el caso ciclico. 



i=l 

+ [vcici + gcic\ + h.c.)] {P + 1) (4.12) 
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El desarrollo hasta aqm es general. De ahora en adelante empezare a 
hacer suposiciones sobre el sistema que voy a estudiar. 



Como se menciono en la seccion 3.0.8 tanto el estado inicial como 



el Hamiltoniano conmutan con P y el estado inicial tiene paridad 
positiva entonces se puede reemplazar P por su valor ya que este 
permanecera invariante. Se tiene entonces: 

n ^ 

^ = 2)-(^4cm + 54cJ+i + /i.c.)] 

i=l 

+2 [vcici + gcic\ + h.c.)] (4.13) 

El termino que no esta en la sumatoria tiene signo opuesto a los 
terminos de interaccion en la sumatoria. En pos de llevar este Ha- 
miltoniano a una forma mas simple, que se pueda escribir como una 
simple sumatoria, hago la primer transformacion lineal q {—lyci 



que cambia los signos de las interacciones. En la figura 14.11 se mues- 
tra esquematicamente como quedan las interacciones de los terminos 
de la sumatoria. Para n par cambian todos los terminos de signo pero 
para n impar el termino que tiene las interacciones entre el el pri- 
mer y el ultimo espin no cambia. En este ultimo caso -n impar- puede 
escribirse entonces el Hamiltonianc0 

n 

H = Y.Hc\c,--) + {vc\c,^, + gc\cl, + h.c.)\ (4.14) 



4.1.2. Transformacion de Fourier Discreta 

Para diagonalizar este Hamiltoniano se realiza primero una trans- 
formada de Fourier discreta de los operadores cj, q dada por 

c] = —^^dle''"^^ ; Wk = 2Tik/n (4.15) 
Se obtiene entonces, para el caso de una cantidad impar de qubits: 



n 

H = 1^(6 + 2v cos Wi)d\di + g sin Wi{d\d\_^ + h.c] 

i=l 



f (4.16) 



^Este Hamiltoniano tambien describe adecuadamente cadenas con n par pero estado inicial im- 
par. 
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• + • 
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I i 
+ + 

• + • 



Figura 4.1: Signo de la integracion entre espines para sistemas de 3, 4, 5 y 6 
espines. En los casos impares queda una interaccion con signo distinto. 

4.1.3. Transformacion de Bogoliubov 

Para terminar la diagonalizacion se propone una transformada 
especial de Bogoliubov de la forma: 



ujaj + Vjan-j 



■- 



dl ■ = Una 



y se lleva a H a la forma diagonal (ver apendice O) 

H = ^[A,(a]a,-^)] + cte 



Xj = \J {h + 2v cos WjY + 4:0^ sin w-j 



(4.17a) 
(4.17b) 

(4.18a) 
(4.18b) 



esto se logra eligiendo v'^) = |[1 ± (6 + 2v cos Wj)/\j\. La constante 
aditiva en el Hamiltoniano no juega ningun papel en el resultado 
final. 



4.1.4. Evolucion Temporal y Valores Medios 

Dado que el Hamiltoniano es diagonal en la base de operadores de 
cuasipartfculas de Bogoliubov es facil calcular la evolucion temporal 
de los valores medios en la representacion de Heisenberg. Esto es 



{0{t)) = Tr{poO{t)) 



dO{t) 
dt 



= i[H,0{t)] 



(4.19) 
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De modo que usando I4.18al y teniendo en cuenta que los ai,aj son 
fermionicos, para aj se tiene 



'dt 



(4.20) 



y analogamente para a!. Integrando elementalmente se tiene 



a]{t) = e'^^'a]{0) y aj{t) = e-'^^'aj{0) (4.21) 

Para poder construir las medidas de entrelazamiento es necesario 
conocer los valores medios de S:^S^_^^ , S:^S^_^^ y S-S-^^. 
Los primeros tres operadores pueden ser escritos como funcion de los 
operadores de cuasiparticulas de Bogoliubov (dejando su dependencia 
temporal impHcita) aplicando las transformadas y considerando que 
como la cadena es ciclica Oj = X^JLi Oi/n. 



n 



t.t 
n-j 



S+S^ 



i+1 



n ^ 



n 



an) 



Q+ C+ 



-E 



n 



a]aj) 



(4.22) 



El valor medio operador restante puede ser escrito en funcion de los 
primeros (Teorema de Wick [|RS41i1 ) como: 



(4.23) 



Para construir los valores medios {S-), {S:^S:^_^-^) y {S:^S:^_^^), como se ve 

en las ecuaciones 1422], es necesario primero conocer {a^jaj), {an-jol^ -} 

y (aja|j_j). Como la dependencia temporal de los operadores de cuasi- 
particulas es simplemente una fase (ec. I4.21D se puede hacer el calcu- 
lo de estos valores medios a t = y luego multiplicar por la fase al 
resumar para armar las expresiones 14.221 Ademas, es crucial notar 
que estos son valores medios de operadores de cuasiparticulas medi- 
dos con respecto al estado inicial, que no es mas que el vacio de los 
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operadores fermionicos de Jordan-Wigner (cj, q). Para poder calcular- 
los es necesario expresar los operadores de Bogoliubov como funcion 
de los de Jordan-Wigner mediante las transformaciones inversas de 
Fourier y de Bogoliubov. Estas son 

aj = Ujd}^ - Vjdn-j ; al_^ = ud\_ - + vdj (4.24) 



d\ = ^y4e-™^-^' (4.25) 



k=i 

De modo que que se obtiene: 



{a]aj) = v] (4.26a) 
(a,_,4_^.) = u] (4.26b) 



a].al_^.) = -u,v,e^^^^'={a^,^,a,Y (4.26c) 



Juntando todo se obtienen los valores medios deseados. 



^ n— 1 < 9 . 2 

^'^^^ " n ^ ^ — ^ ~ ^ (4.27a) 
{StS-,^ = 1 g ^g^^'"^ (4.27b) 

TL A- 



ism2Xjt] (4.27c) 

(^■^m) = (^•)' + |(^.^^ii)|'-|(^r^r+i)r (4.27d) 



4.2. Entrelazamiento 
4.2.1. Generalidades 

En este caso la funcion a{t) que determina el entrelazamiento de 



uno con el resto es simplemente {S-) + 1/2 (ver ecuacion 1427aD. Ana 



lizaremos primero la forma de esta y daremos resultados generales 
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y luego indicaremos generalidades sobre el entrelazamiento de pares 
que es una funcion mucho mas complicada de la ecuaciones 14.271 
La envolvente superior de a{t) se obtiene haciendc@ sin^ A^^ = 1 V j y 
es 

■1 ,2^;_2 



1 



= -E 



4g sin w. 



{b + 2v cos Wj^ + 45f2 sin Wj 



(4.28) 



Este resultado es una generalizacion natural de los anteriores. Es 
una suma de gausianas con centros en 6 = —2vcoswj y de ancho 
25^1 sin Wj\. Para n par la cantidad de picos distintos es n/2+1, mientras 
que para n impar, lo que nos interesa ahora, hay (n — 1) /2 picos. Como 
siempre, nos interesa saber cuando esta funcion es menor o mayor a 
1/2 ya que esto determina si se llega al entrelazamiento maximo o 
no. 

La condicion para que un pico j aislado llegue a un valor mayor 

a 1/2 es (6 + 2v cos w-,)^ < ig'^ sin^ Wj de lo que se puede deducir que si 



V ^ gyb <2v se tiene a{t) > 1/2. Estos resultados pueden constatarse 



comparando las figuras [4.2(a)| y [4y3(a)| con las |4.2(b)1 y |4.3(b)| donde se 
muestra el entrelazamiento de uno con el resto para g = v = 1 (hay 
saturacion) y g = v/2 = 1 (no hay saturacion) correspondientemente, 
para dos sistemas de distinta cantidad de qubits. 

Si g <vno solo se deja de llegar a la saturacion sino que a medida 
que g se hace mas chico que v se empiezan a distinguir los picos. La 
distancia entre dos picos consecutivos es 6b = 2v \ coswj 



2v\^^^^p-\6wj = ^1 sin Wj \ y el ancho de un pico es 2g\ sin Wj\. Podemos 
estimar que se veran picos claros cuando la separacion entre ellos sea 
mayor que sus anchos; eso es 5^ < 27rv/n. Esto puede constatarse con 
los graficos [4.21 y 14.31 Para el caso de n = 5 se tiene 2TTv/n « 1,26 por 
lo que en todos los casos mostrados se deberlan distinguir picos sin 
embargo no es asi porque para v = 1 la. saturacion borra los picos. 
Para n = 15 se tiene 27rv/n « 0,42 y como se ve solo se resuelven 



cos Wj+i\ 



los picos en los graficos |4.3(c)| y |4.3(d)| que corresponden a. v = 0,3 y 
V = 0,1 respectivamente; en |4.3(a)| y |4.3(b)| que corresponden a. v = 1 
y V = 0,5 no se resuelven picos. 

En el llmite asintotico de 6 > {v,g) se obtiene a„^ « 2g^/b'^ que es 
el mismo resultado que para tres qubits y por ende para cualquier n 
impar, en este Ifmite, Ci « 2\/2g/\b\. 



^Solo es posible que se de esta situacion si las razones entre las frecuencias no son niimeros 
racionales. 
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(c) g = 0,3 (d) g = 0,1 



Figura 4.2: Entrelazamiento maximo de uno con el resto Ci y de a pares C2 para 
= 5 en funcion de 6/f para distintos valores de g (en unidades de v). La Imea 
gruesa es el entrelazamiento de uno con el resto Ci. Lo que se ve por debajo de esta 
es el entrelazamiento de pares C2 para muchos tiempos distintos. 



En general el entrelazamiento de pares esta muy por debajo del 
entrelazamiento de uno con el resto y lejos de la resonancia es practi- 
camente una funcion monotonamente decreciente del campo 6; en el 
limite asintotico de 6 > {v.g) se tiene C2 ~ '^g/h. Sin embargo, para 
los casos en que hay resonancia la concurrencia tambien resuena y 
veremos mas adelante que se pueden encontrar resultados interesan- 
tes. 



4.2.2. Entrelazamiento y Deeoherencia 

La evolucion temporal para sistemas de muchos espines es en ge- 
neral bastante complicada, pero hay casos para los que se pueden 
encontrar resultados relevantes, en los otros indicaremos generali- 
dades. 



Empezamos con los complicados. La figura 14.41 muestra la evolu 



cion temporal del entrelazamiento en cadenas de 5 y 15 qubits pa- 
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(c) g = 0,3 (d) g = 0,1 



Figura 4.3: Entrelazamiento maximo de uno con el resto Ci y de a pares C2 para 
n = 15 en funcion de b/v para distintos valores de g. La Imea gruesa es el entrelaza- 
miento de uno con el resto Ci. Lo que se ve por debajo de esta es el entrelazamiento 
de pares C2 para muchos tiempos distintos. 



ra un campo magnetico con el que se llega a la saturacion h = 0,5; 
V = g = 1 (ver figs |4.2(al1 y |4.3(a)| ). Se observa que para 5 qubits la de- 
pendencia temporal del entrelazamiento de uno con el resto es cuasi- 
periodica y recorre continuamente todos los valores desde cero hasta 
la saturacion mientras que para 15 qubits esta periodicidad parece 
desaparecer y el sistema no recorre todos los valores posibles del en- 
trelazamiento de uno con el resto sino que oscila entre la saturacion 
y una cota minima. En realidad, ambos comportamientos son cua- 
siperiodicos y a un tiempo suficientemente largo es natural esperar 
que, para cualquiera de los casos, se de un momento en que el entre- 
lazamiento de uno con el resto se anule. Sin embargo, para n = 15 el 
entrelazamiento de uno con el resto permancece principalmente por 
arriba de una cota. El entrelazamiento de pares, en cambio, es menor 
y ocurre para tiempos mas esporadicos cuando crece la cantidad de 
qubits en la cadena. 

El entrelazamiento de pares, de algun modo, es un indicador de 
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u 




u 



t-vjn 
(a) n=5 




Figura 4.4: Entrelazamiento de uno con el resto Ci (linea roja solida) y de a pares 
C2 (linea azul a trazos) en funcion del tiempo en unidades de h/v con b = 0,5, g = I 
(saturan Ci y C2) para n = 5 y n = 15. 



que "tan cuantico" es un subsistema. Al agrandarse la cadena la can- 
tidad de espines que esta monitoreando □ a los dos espines estudia- 
dos se ve que la probabilidad de que se genere entrelazamiento en el 
subsistema es menor, es decir el sistema es "menos cuantico". Resu- 
miendo, y a groso modo, se puede decir que, al agrandar el tamano 
del sistema, un subsistema de un tamano dado pierde la capacidad de 
que sus partes esten entrelazadas entre si, pero estas permancecen 
constantemente entrelazadas con el resto del sistemaH . Esta obser- 
vacion es una de las ideas que dio a lugar a la teorfa de la decohe- 
rencia que da un gran paso en la explicacion de como, a pesar de que 
el mundo responde a las leyes e la mecanica cuantica, lo percibimos 
como clasico (ver apendice[D]). 



^ Inter actuando con, haciendoles de bano termico u observando constantemente. 
*E1 resto del sistema viene a estar haciendo de medio ambiente o bafio termico para el subsitema 
que se esta estudiando. 
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Figura 4.5: Entrelazamiento de uno con el resto Ci (Imea roja solida) y de a pares 
C2 (linea azul a trazos) en funcion del tiempo en unidades de h/v con b = A, g = 1 
(no satura Ci ni C2) para n = 5 y n = 15. 



La figura 14.51 muestra la evolucion temporal del entrelazamiento 
para sistemas con la misma cantidad que qubits que antes pero con 
un campo magnetico mas fuerte que hace que el comportamiento es- 
te fuera de la saturacion b = 4; v = g = 1 (ver figs |4.2(a)| y |4.3(a)D . 
El comportamiento, para el entrelazamieno de uno con el resto, es 
similar al anterior, pero en este caso el maximo alcanzado no llega a 
saturar. El entrelazamiento de pares tambien disminuye al agrandar 
la cadena, sin embargo, es mayoij^ que en el caso anterior. El incre- 
mento en campo magnetico disminuye el entrelazamiento global (de 
uno con el resto) impidiendo que las interacciones locales tengan in- 
fuencia a largo rango. Contrariamente se hacen mas significativas 
las caracteristicas que solo dependen del corto rango. Para un campo 
magnetico suficientemente grande el entrelazamiento de pares pacti- 



^Es mayor en el sentido de que, en promedio, hay mas tiempos en el que el entrelazamiento de 
pares esta per encima de una cota. 
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camente "sigue" al entrelazamiento de uno con el restcH del mismo 
modo que se habia observado para el caso de solo tres espines (ver 
figura[33(d)|). 



4.2.3. Resonancias 



En el caso en que hay resonancia, ver figuras |4.2(d)| y |4.3(d)l la 
dependencia temporal es muy simple. La resonancia se da cuando la 
anisotropia de la interaccion es muy chica g <^ v. Los campos para 
los cuales el sistema resuena cumplen b = —2v cos Wj de modo que los 
unicos terminos relevantes en a{t) son justamente el j y el n — j, que 
corresponden ambos al mismo pico. En tal caso se puede aproximar 
a{t) como: 

2 

a{t) = - siia^ Xjt ; Xj = 2g\smwj\ (4.29) 



La figura |4.6(a)| muestra la evolucion temporal del entrelazamiento 
para un sistema de 15 qubits con v = 1 y g = 0,1. Se muestran re- 
sultados para el campo b justo en la resonancia del pico de j = 7, a 
1/10 de la distancia entre este pico y el anterior y justo en el medio de 
estos dos picos 0. Es notable la sensibilidad a pequenas variaciones 
de campo. Por mas que en este caso el entrelazamiento no llegue a 
saturar, a los fines practicos de construir una computadora cuantica 
este caso es de gran importancia ya que no solo la dependencia tem- 
poral es simple sino que con pequenas variaciones de campo puede 
"prenderse" o "apagarse" facilmente la formacion de entrelazamien- 
to. Tambien es importante notar que el unico autovalor relevante del 
sistema Xj = 2g\ sin Wj\ se hace muy pequeno de modo que el tiempo 
necesario para llegar al maximo se hace muy grande en comparacion 
con los tiempos naturales del sistema fuera de la resonancia. 

Cerca de la resonancia la concurrencia de pares para estados con 
paridad negativa tampoco es dificil de calcular. En particular se tiene 
que en este limite las ecuaciones 14.271 pueden aproximarse como: 



2 

(5+ 5,^1) « -| cos Wj\ sin^ Xjt (4.30) 
(5+5'+ 1) « -i|sinw;^ sin2A^i| (4.31) 

T7/ 



^En estos casos de muchos qubits, sin embargo, aparecen comportamientos mas complicados 
como saturaciones y "revivals". 

^El pico j = 7 es el pico de mas a la derecha en la figura |4r3(a)t el anterior j = 6 es el segundo de 
derecha a izquierda. 
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Figura 4.6: Entrelazamientos Ci y C2 en funcion del tiempo en unidades de h/v 
en resonancia para el pico j = 7 de un sistema de 15 qubits con = 0,1. (a) Ci para 
campo resonante, a 1/10 de la distancia entre picos y justo entre picos. (b) C2 para 
campo resonante y a 1/10 de la distancia entre picos (ver texto). 



de modo que los picos aparecen para algunos de los campos resonante 
en los tiempos en que sin Xd = 1 en los cuales la altura de estos es 



Co = - 



cos w 



n 



Ismw 



2 2 

— V r COS^ Wj 

n n'^ 



(4.32) 



Estos picos son de un valor mucho menor al del entrelazamiento de 
uno con el resto y aparecen adyacentes a los tiempos en los que este 
se anula y son de paridad opuesta al estado inicial. Se muestran en 
la figura [476(b)| los resultados para el campo 6 justo en la resonancia, 
a 1/10 de la distancia entre picos se ve que la sensibilidad del entre- 
lazamiento de pares frente a pequenas desviaciones de la resonancia 
es aun mayor que para el entrelazamiento de uno con el resto. 
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4.2.4. Tiempos cortos 



Por ultimo analizo el extremo v = 0. En este caso todos lo picos 
colapsan en el origen. Para 6 = se tiene ademas que a™ = (n — l)/n 
de modo que hay un tiempo para el cual se llega a la saturacion. 
Lo interesante de este caso en particular es que es facil hallar este 
valor de tiempo. El procedimiento consiste simplemente en hacer un 
desarrollo en serie de a{t) de un orden suficientemente grande tal 
que siga a la funcion hasta que esta vale 1/2. El orden necesario del 
desarrollo para que cumpla esta condicion es 10 (como puede verse 
en la fig 14.71) y el desarrollo que se obtiene, curiosamente, no depende 
de n. 

ait) « 2gH' - 2gH' + -gH' - -gH' + —g^'t^' (4.33) 
El tiempo en el que se alcanza la saturacion a{t) = 1/2 es t « 0,6012/^'. 





(a) n 



(b) n = 15 



Figura 4.7: Entrelazamientos de uno con el resto Ci en funcion del tiempo en 
unidades de h/g y serie de Taylor en el caso que b = v = para n = 5 y n = 15. 
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Conclusion 



Se examino la formacion de entrelazamiento a partir de un estado 
completamente alineado en una cadena XY con un campo magneti- 
co uniforme. Contrariamente a lo que se podna haber esperado, el 
comportamiento no es monotono con el campo magnetico. El maximo 
entrelazamiento alcanzado presenta resonancias y "plateaus". 

Los "plateaus" ocurren para anisotropias que son del mismo orden 
mayor que el parametro de "hopping" y saturan cuando es mayor. Se 
mostro, ademas, que un aumento de la anisotropia no necesariamen- 
te garantiza una prolongacion del tiempo durante el cual el sistema 
permanece maximamente entrelazado. 

Las resonancias ocurren cuando la anisotropia es pequena en com- 
paracion con el parametro de "hopping". Son muy marcadas para el 
entrelazamiento de uno con el resto y mas aun para el entrelazamien- 
to de pares. Se encontro, ademas, que estas tienen un valor finito para 
anisotropias arbitrariamente chicas pero no nulas. Las resonancias 
permiten identificar configuraciones aptas para generar o no generar 
entrelazamiento en circuitos cuanticos con pequenas variaciones del 
campo magnetico. 

Estos resultados y otros han sido aceptados para su publicacion 
en la revista Physical Review A en un trabajo conjunto con Raul Ros- 
signoli. Entre los resultados que se incluyen figuran ademas: cadenas 
con cantidad par e impar de qubits, resonancia de pares entrelazados 
de paridad positiva, un estudio detallado del comportamiento para 
V = g y b = donde la evolucion es estrictamente periodica y resulta- 
dos de diagonalizacion directa para n = 4. 
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Apendices 
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Apendice A 

La descomposicion de 
Schmidt 



El teorema de Schmidt es una herramienta muy poderosa para el 
analisis de sistemas bipartitos. Explico el teorema sin demostraci6n0, 
menciono algunos corolarios inmediatos y hago un comentario. 

Teorema 

Sea un estado de un sistema bipartite de dimension arbitraria 
descripto per un ket Existen estados ortonormales |u) para el 
sistema A y estados ortonormales para el sistema B tales que: 



donde los son numeros reales no negativos que satisfacen ^ • = 1 
conocidos como coeficientes de Schmidt. La cantidad de coeficientes de 
Schmidt no nulos se denomina numero de Schmidt. 

Corolarios 

Algunos corolarios que pueden demostrarse inmediatamente son: 

■ Los autovalores de las matrices reducidas de cada subsistema son 
iguales. 

■ Si un estado es separable su numero de Schmidt es 1. 

^La demostracion se basa en la Descomposicion en Valores Singulares (DVS) y puede verse en el 
libro de Nielsen y Chuang INCOOJ . 




(Al) 
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■ El numero de Schmidt no cambia si se hacen operaciones unita- 
rias sobre alguno de los subsistemas. 

Comentario 

A pesar de su potencia para revolver problemas la descomposi- 
cion de Schmidt vale solo para un par de subsistemas de dimension 
arbitraria. Es decir, cada subsistema puede estar compuesto por mas 
de un elemento pero no pueden siempre realizarse descomposiciones 
del tipo de Schmidt. Por ejemplo, para un sistema de tres componen- 
tes A,B,C pueden realizarse descomposiciones del tipo Xi Iiab) l^c), 
J2i I^a) liBc) K liAc) lis) donde los kets de cada subespacio son 
ortonormales y los numeros reales positives. En cambio no necesa- 
riamente pueden realizarse descomposiciones del tipo {ib) \ic) 

que cumplan que los kets de cada subespacio sean ortonormales y los 
Xi sean numeros reales no negatives. 



Apendice B 

Paridad 



B.l. Paridad y Subsistemas 

Demostrare que si un estado tiene paridad definida, la matriz den- 
sidad reducida de cualquier subsistema conmuta con el operador pa- 
ridad. P \iIjab) = ±\i^)^ [P: TrBPAB] = 0. 



Definiciones 

Defino el operador paridad general P con autovectores |+),|— ) de 
manera tal que P |+) = |+) y P |— ) = — |— ). 
Un estado tiene paridad definida si es autoestado de P. 



Demostracion 

En general podra haber varies vectores — «)) en una base que 
sean autovectores de P. En particular, si estos vectores forman una 
base completa del espacio del Hilbert, cualquier estado podra ser es- 
crito como: 



Para un sistema bipartite esto se extiende a: 

\^ab) = ^a,bj 1+,) |+_^) + ^g^hj |-,) 

ij ij 

+ J2c^dj 1+,) I-,) + ^ej, I-,) 1+,) (B.2) 

ij ij 
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Capitulo B. Paridad 



Es claro que los dos primeros terminos tienen paridad positiva y los 
restantes negativa. Constmyo entonces dos autovectores con paridad 
definida y vere si al tomar la traza parcial sobre un subsistema se 
conserva la paridad. Estos son (juntado las sumas): 

I^Jb) = X] ^^'^^ + Skhi \-k-i)} 

ijkl 

I^ab) = + ^'^•^^ \-k+i)} (B.3) 

ijkl 

La matriz densidad para paridad positiva = l^'J^) {^ab\ 
Pab = X] {^^^AK + 9khiglhp \-k-i) {-o-p\ 

ijklmnop 

+aibjglhl {-m-n\ + QkhalXi {+m+n\] (B.4) 

Al tomar la traza parcial sobre un subsistema {B) se obtiene la matriz 
reducida: 

Pa= X {^^^j^mK (+m| + Qkhglhl \-k) (o-|} (B.5) 

ijklmnop 

Se ve claramente que [p\, P] = 0. La demostracion para paridad ne- 
gativa es exactamente igual. QED. 

En realidad podria haberse demostrado lo anterior simplemente di- 
ciendo que como la traza parcial implica una reduccion que deja so- 
lamente elementos que son diagonales en las partes trazeadas si la 
matriz densidad del sistema conmuta con P tambien tiene que hacer- 
lo la matriz reducida. Estos resultados tambien son validos si pab es 
un estado mixto. 



Caso n = 1; un subsistema de dos niveles de distinta paridad 

En este caso, por ser 2 la dimension del subsistema que se desea 
estudiar, la matriz reducida es la ecuacion. IB.5I pero sin las sumas. 



Pa 



+ + + 



bien en forma matricial (teniendo en cuenta que Trpi = 1): 



Pi 



a 
l-a 



(B.6) 
(B.7) 



(B.8) 



B.2. Paridad de espm 
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Caso n = 2; dos subsistemas de dos niveles de distinta paridad 

En este caso la dimension del subespacio es 4. Los autoestados de 
P son: 





- I++>,I 


1-.) - 





de modo que la matriz densidad sera (directamente en forma matri- 
cial): 



/ 


a 








d* \ 







b 


e* 










e 


c 





v 


d 








// 



(B.IO) 



con f = l — a — b — c debido a que Trp2 = 1. 



B.2. Paridad de espin 

En un sistema de espines 1 /2 puede definirse un operador -paridad 
de espiiA que es un buen operador paridad asociado a la orientacion 
de los espines. 

p = Sf+n/2) ^ -Q gZ^(5f +n/2) (B.ll) 

i 

donde n es la cantidad de espines totales en el sistema. 
Asi definido este sera un buen operador paridad definido sobre la ba- 
se de 5^ Esto es la base |T),U) donde |T) = | IT),'?" \i) = -\ U). 
Se ve que este es efectivamente un buen operador paridad estudian- 
do como actua sobre algunos estados. En particular sobre el estado 
||^«-T ii^^-i es inmeditado ver que: 

p ||^®"T _ (^_]_^'^T ||^®"T I (B.12) 

Puesto en palabras: el operador paridad distingue si hay una canti- 
dad par impar de espines para "arriba". 



En ingles se los suele llamar phase flip symmetry. 
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Apendice C 



Transformada especial 
de Bogoliubov ■ BCS 



En este apendice muestro como diagonalizar un Hamiltoniano es- 
pecifico utilizando una transformada especial de Bogoliubov, tambien 
llamada tranfsormada de BCS por Bardeen, Cooper y Shrieffer, que 
la usaron en su famoso paper de 1957 en el cual presentaron la pri- 
mer teoria puramente cuantica de la superconductividad. 

Sea un sistema de fermiones cuyos operadores de creacion y aniqui- 
lacion son di y dj descripto por un Hamiltoniano del tipo 



n 

H = ^ [x,, {d\di + d\_^dn-^ + yi (d\d\_^ + dn-id^ 
1=1 



(C.I) 



Un transformacion especial de Bogoliubov de operadores de cuasi- 
particulas fermionicas y a] viene dada por 



d] = u^a] + v^ttn-i 
dl_i = Uial_^ - Vitti (C.2) 



donde se debe tener que Ui y Vi reales y uf + vf = 1 para que los opera- 
dores de cuasiparticulas sean fermionicos. Haciendo los reemplazos 
correspondientes y agrupando terminos con iguales operadores se ob- 
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tiene: 

n 

H = ^a-fli [xi (uj - Vi) - 2yiUiVi] 

i=l 

+al_^an-z [xi [u] - vl) - 2y,u^v^] 

+a|j_^aj [2xiUiVi + yi [u] - vl)] 

+aian-i [2xiUiVi + yi (u- - v^)] 

+ - yiUiVi + XiVi (C.3) 

Como se desea lleva el Hamiltoniano a una forma diagonal se debe 
exigir que 2xiUiVi + y,,{uf — vf ) se anule. Esta ecuacion junto con la 
exigencia de que los operadores de cuasipartfculas sean fermiones 
(u| + v1 = 1) tiene como posibles soluciones 

de modo que el Hamiltoniano se escribe como: 

H = [a\a^ + a\_^an-^ + ^ (C.5) 



Apendice D 

Decoherencia 



Con el fin de dar alguna idea sobre la decoherencia incluyo en 
este apendice una traduccion que hice de la introduccion a un curso 
[IPZ00I1 que dieron Paz y Zurek sobre Decoherencia inducida por el 
medio ambiente y la transicion cuantico-clasica. 

El origen cuantico del mundo clasico ha sido tan diffcil de imagi- 
nar por los padres de la mecanica cuantica que llegaron a postular 
su existencia independiente, como lo hizo Niels Bohr, o trataron de 
buscar explicaciones mas complejas con fundamentos clasicos, como 
lo hicieron de Broglie y Einstein. El problema radica en el principio 
cuantico de superposicion que, en efecto, expande el posible espectro 
de estados posibles a una infinita superposicion de ellos. Por ende 
superposiciones coherentes de gatos muertos y vivos tienen - dentro 
del formalismo cuantico - el mismo derecho de existir que cualquiera 
de sus dos alternativas clasicas (vivo o muerto). Dentro de todos los 
estados cuanticos posibles, los clasicamente permitidos son excepcio- 
nales. La cantidad de estados posibles para un sistema cuantico es 
enorme comparado con la cantidad que uno encuentra en sistemas 
clasicos. Sin embargo es un hecho veridico que solo encontramos ob- 
jetos clasicos en un subgrupo muy pequeno de todos los estados po- 
sibles (en principio permitidos). Debe uno explicar el origen de esta 
aparente regla de "super-seleccion" que previene la existencia de la 
mayorfa de los estados posibles de un dado sistema. La decoherencia 
y su principal consecuencia, la super-seleccion inducida por el medio 
ambiente -einseleccion-, dan cuenta de este hecho experimental del 
mundo fisico. 

La causa de la decoherencia es la interaccion entre el sistema y 
su medio ambiente. Bajo una variedad de condiciones, que son muy 
faciles de satisfacer para un objeto macroscopico, se da la einselec- 
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cion, de un pequeno subgrupo de estados cuasi-clasicos dentro de la 
enorme cantidad de estados cuanticos posibles. La clasicalidad es, 
entonces, una propiedad emergente inducida por el sistema y su in- 
teraccion con el medio ambiente. Las superposiciones arbitrarias de 
estados permi tides son dejadas de lado en boga de un subgrupo prefe- 
rido de "estados punteros" que emerge. Estos estados preferidos son 
los candidates a los estados clasicos. 

El rol de los procesos de decoherencia en inducir la clasicalidad 
se ha aclarado solo recien en la ultimas tres decadas. La idea cen- 
tral es simple: El medio ambiente de un sistema cuantico puede, en 
efecto, monitorear sus estados mediante una continua inter accion. La 
marca que los estados dejen en el medio ambiente contendra la infor- 
macion de cuales seran los elegidos. Los estados que prevalecen son 
los que puedan interactuar con el medio ambiente sin ser perturba- 
dos. La propiedad clave que estos estados cuasi-clasicos tendran es, 
entonces, su insensitividad al monitoreo y consecuentemente su re- 
sistencia al entrelazamiento causado por la interaccion con el medio 
ambiente. Estos son tambien, casi por definicion, los unices estados 
que describen al sistema por si solo y no conjuntamente con el medio 
ambiente. Los otros estados son superposiciones entre el sistema y 
medio ambiente y solo conservan su pureza cuando ambos son anali- 
zados conjuntamente en un mayor "super-sistema". 

El hecho que la interaccion entre sistemas cuanticos produce en- 
trelazamiento es bien conocido desde los principios de la mecanica 
cuantica. En efecto como las ideas de la decoherencia y la einselec- 
cion se basan en la teoria cuantica, y solamente en la teoria cuantica, 
es licito preguntarnos porque se tardo tanto en llegar a una expli- 
cacion natural de los origenes cuanticos de la clasicalidad. Hay va- 
rias posibles explicaciones para esta demora. El prejuicio que parece 
haber demorado mas el estudio de sistemas cuanticos abiertos (que 
iteractuan con un medio ambiente) esta basado en un mode clasico 
de pensar sobre el Universe. Dentro del contexto de la fisica clasica 
todas las preguntas siempre se han podido responder pensando en 
un sistema cerrado (sin interaccion con el medio ambiente). La estra- 
tegia basica consiste en asegurar la aislacion del sistema involucra- 
do agrandando el sistema a estudiar (incluyendo su medio ambiente 
mas proximo). La expectativa era que de este mode uno podia redu- 
cir cualquier sistema abierto a uno mas grande pero cerrado. Esta 
estrategia funciona en la fisica clasica donde el entrelazamiento pue- 
de ayudar a satisfacer leyes de conservacion como la de la energia. 
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Falla en el caso cuantico bajo consideracion ya que en este caso es la 
informacion la que se desea que no se esparza. Es mucho mas dificil 
contener informacion cuando un sistema se vuelve mas grande. En 
fin, el unico sistema macroscopico verdaderamente aislado es el Uni- 
verso. Y nosotros, los observadores, ciertamente no estamos en con- 
diciones de estudiarlo desde afuera. 
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